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ABSTRACT 

Hypervelocity stars (HVSs) are stars ejected completely out of the Milky Way by three-body in- 
teractions with the massive black hole in the Galactic center. We describe 643 new spectroscopic 
observations from our targeted survey for HVSs. We find a significant (3.5a) excess of B-type stars 
with large velocities +275 < v r f < 450 km s _1 and distances d > 10 kpc that are most plausibly 
explained as a new class of HVSs: stars ejected from the Galactic center on bound orbits. If a Galactic 
center ejection origin is correct, the distribution of HVSs on the sky should be anisotropic for a survey 
complete to a fixed limiting apparent magnitude. The unbound HVSs in our survey have a marginally 
anisotropic distribution on the sky, consistent with the Galactic center ejection picture. 
Subject headings: Galaxy: halo — Galaxy: center — Galaxy: stellar content — Galaxy: kinematics 
and dynamics — stars: early-type 



1. INTRODUCTION 

HVSs are a natural consequence of a massive black 
hole (MBH) in a dense stellar environment like that in 
the Galactic center. Hills (1988) first pointed out that a 
stellar binary encountering the Milky Way's central MBH 
can eject one member of the binary as a HVS traveling 
at ^1,000 km s _1 . HVSs differ from classical "runaway" 
stars because 1) HVSs are unbound and 2) the classical 
supernova ejection (Blaauw 1961) and dynamical ejection 
(Poveda et al. 1967) mechanisms that explain runaway 
stars cannot produce ejection velocities larger than 200 - 
300 km s^ 1 for B-type stars (Leonard 1991, 1993; Porte- 
gies Zwart 2000; Davies et al. 2002; Gualandris et al. 
2004; Dray et al. 2005). The first HVS discovered, by 
comparison, has a heliocentric radial velocity of +853 
km s _1 and a Galactic rest-frame velocity of at least 
+709 + 12 km s" 1 (Brown et al. 2005a), many times that 
needed to escape the Milky Way. Photometric follow-up 
shows that the object is a slowly pulsating B main se- 
quence star (Fuentes et al. 2006). Only interaction with 
a MBH can plausibly accelerate this 3 M Q main sequence 
B star to such an extreme velocity and distance (^110 
kpc). 

HVSs are fascinating because they can be used to un- 
derstand the nature and environs of MBHs (see the re- 
cent theoretical work of Gualandris et al. 2005; Holley- 
Bockelmann et al. 2005; Ginsburg & Loeb 2006a; Levin 
2006; Perets et al. 2006; Baumgardt et al. 2006; Sesana 
et al. 2006; Bromley et al. 2006; Merritt 2006; Demarque 
& Virani 2007; Ginsburg & Loeb 2006b; O'Leary & Loeb 
2006; Gualandris & Portegies Zwart 2007; Kollmeier 
& Gould 2007). The trajectories of HVSs also pro- 
vide unique probes of the shape and orientation of the 
Galaxy's dark matter halo (Gnedin et al. 2005). Discov- 
eries of additional HVSs (Edelmann et al. 2005; Hirsch 
et al. 2005; Brown et al. 2006a, b) are starting to provide 
suggestive limits on the stellar mass function of HVSs, 
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the origin of massive stars in the Galactic Center, and 
the history of stellar interactions with the MBH. Clearly, 
a larger sample of HVSs will be a rich source for further 
progress on these issues. 

We have designed a successful targeted survey for new 
HVSs. We use the 6.5m MMT and the Whipple 1.5m 
Tillinghast telescopes to obtain radial velocities of faint 
B-type stars, stars with lifetimes consistent with travel 
times from the Galactic center but which should not oth- 
erwise exist in the distant halo. Four earlier HVS discov- 
eries from this survey are published elsewhere (Brown 
et al. 2006a, b). Here, we present spectroscopic observa- 
tions of 643 new HVS candidates. 

Our paper is organized as follows. In §2 we discuss our 
target selection and spectroscopic identifications. In §3 
we present evidence for stars ejected from the Galactic 
center on bound orbits. In §4 we show that the unbound 
HVSs have a marginally anisotropic distribution on the 
sky. We conclude in §5. Our new observations are listed 
in Appendix A. 

2. DATA 
2.1. Target Selection 

Brown et al. (2006b) describe our target selection for 
candidate HVSs. Briefly, we use Sloan Digital Sky Sur- 
vey (SDSS) photometry to select late B-type stars by 
color. Figure 1 presents the HVS survey color-color se- 
lection. We use the 6.5m MMT telescope to observe ev- 
ery star with 17 < g' < 19.5 in this color range. We have 
observed 136 new HVS candidates and have completed 
the high declination region of the SDSS Data Release 4 
(DR4, Adelman-McCarthy et al. 2006). The MMT sur- 
vey is now 80% complete across SDSS DR4 and covers 
an area of ^5,000 deg 2 or 12% of the entire sky. 

We have also carried out a new, complementary HVS 
survey targeting bright stars with 15 < g' Q < 17 and 
B-type colors in the SDSS. We use the 1.5m Tilling- 
hast telescope and the FAST spectrograph (Fabricant 
et al. 1998) to observe these bright HVS candidates. 
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Fig. 1. — Color-color diagram showing the target selection re- 
gions for the MMT (dashed line) and FAST (solid line) samples, 
and the spectroscopic identifications for the combined set of 894 
HVS candidates. The brighter FAST sample (15 < g' Q < 17) tar- 
gets bluer stars and finds a number of stars with B3 and B4 spectral 
types (solid dots). 



We are able to survey efficiently much bluer stars with 
FAST because contamination from white dwarfs is less 
problematic at bright magnitudes. Figure 1 shows the 
FAST color-color selection: 0.1 < (u' - g') < 1.5 and 
(10.67( 5 '-r') +3.5) < {u'-g') < (10.67(.g'-r') +5.07). 
This region follows the stellar sequence of B-typc stars 
in the SDSS photometric system (Fukugita et al. 1996). 
Following Brown et al. (2006a), we impose a color cut 
—0.5 < (r' — i')o < to reject objects with non- 
stellar colors. We also exclude the region of sky between 
b < —1/5 + 50° and b > 1/5 — 50° to avoid excessive 
contamination from Galactic bulge stars. There are 746 
SDSS DR4 candidate B stars with 15 < g' < 17 in the 
selection region. We have observed 514 stars, 69% of 
this total. The average surface number density of bright 
HVS candidates is 0.13 deg~ 2 . Thus we have surveyed 
an effective area of ~4000 deg 2 for HVS candidates with 
15 < g' a < 17 with FAST. 

2.2. Spectroscopic Observations and Radial Velocities 

New observations at the 6.5m MMT telescope were ob- 
tained with the Blue Channel spectrograph on the nights 
of 2006 May 23-26 and 2006 June 19-20. The spectro- 
graph was operated with the 832 line mm" 1 grating in 
second order and a 1.25" slit. These settings provided 
a wavelength coverage of 3650 A to 4500 A and a spec- 
tral resolution of 1.2 A. Exposure times ranged from 5 
to 30 minutes and were chosen to yield S/N = 15 in the 
continuum at 4000 A. Comparison lamp exposures were 
obtained after every exposure. 

Observations at the 1.5m Tillinghast telescope were 
obtained with the FAST spectrograph over the course of 
19 nights between 2006 January 1 and 2006 July 21. The 
spectrograph was operated with a 600 line mm -1 grating 
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Fig. 2. — Spectrum of the z = 2.43 quasar located at 
15 h 58 m 5rf85, +22 d 21 m 59!9 (J2000). 



and a 2" slit. These settings provided a wavelength cov- 
erage of 3500 A to 5400 A and a spectral resolution of 2.3 
A. Like the MMT observations, exposure times ranged 
from 5 to 30 minutes and were chosen to yield S/N =15 
in the continuum at 4000 A. Comparison lamp exposures 
were obtained after every exposure. 

Spectra were extracted using IRAF 1 in the stan- 
dard way. Radial velocities were measured using 
the cross-correlation package RVSAO (Kurtz & Mink 
1998). Brown et al. (2003) describe in detail the cross- 
correlation templates we use. The average uncertainty is 
±11 km s" 1 for the B-type stars. 

2.3. Selection Efficiency and Completeness 

Our combined observations total 894 HVS candidates 
with 15 < g' < 19.5, of which 699 (78%) are stars of 
B spectral type. The remaining objects include 185 DA 
white dwarfs, 3 other (DZ, DB) white dwarfs, 1 com- 
pact blue galaxy (see Kewley et al. 2006), 1 quasar, and 
5 objects of uncertain identification, possibly broad ab- 
sorption line quasars. The one obvious quasar is a previ- 
ously unidentified g' = 19.67 object at z — 2.43, located 
at 15 h 58 m 51!85, +22 d 21 m 59!9 (J2000). Figure 2 presents 
its spectrum. The quasar has colors (u' — g') = +0.80 
and (g' — r') = —0.28, unusually red in (u 1 — g') and 
blue in (</' — r') compared to other quasars at similar 
rcdshifts (Schneider et al. 2005). 

Figure 1 plots the spectroscopic identifications of the 
HVSs, other B-type stars, and the white dwarfs. Inter- 
estingly, we find stars with spectral types as early as B3 
in the bright FAST sample: the solid dots in Figure 1 are 
the 61 stars with spectral types B3 and B4. If these early 
B-type stars are main sequence stars, they are located at 
surprisingly large heliocentric distances 30 < d < 60 kpc. 
We discuss the early B-type stars further in §3. 

Our survey is now 100% complete across the high 
declination-region of SDSS DR4 located at 25 < I < 
210°, b > 0° in Galactic coordinates. Curiously, we find 
no new unbound HVSs in the completed region. Our 
previous discovery rate with the MMT was 1 HVS per 
~50 B-type stars (Brown et al. 2006b), yet our new set 
of observations contains 245 (57 from the MMT) B-type 
stars in the complete region. The absence of an unbound 
HVS in this large region suggests that the distribution 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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of HVSs on the sky may have structure. We address the 
spatial distribution of HVSs in §4. 

For the remainder of this paper, we consider a "clean" 
sample containing 560 stars selected only from the over- 
lapping region of color-color space: —0.38 < (g' — r') < 
-0.28 and 2.67(. 9 ' - r')o + 1.30 < (u' - g') < 2.67($' - 
r')o + 2.0. Main sequence B stars in this selection re- 
gion have masses of 3 - 5 M© . We also consider the full 
sample of 699 B-type stars at all colors. The full sample 
contains earlier B-type objects from the FAST sample 
possibly including main sequence stars with masses up 
to 7 Mq. We emphasize that our photometric selection 
is complete and can detect a B-type star at any velocity. 

3. BOUND HYPERVELOCITY STARS 

The observed distribution of distances and velocities 
of our B-type stars provides ample evidence that they 
are not main sequence run-aways. Curiously, the veloc- 
ity distribution reveals an asymmetry of stars with large 
positive radial velocities. Because the large positive ve- 
locity outliers are very unlikely main sequence run-aways, 
we argue that the stars represent a class of bound HVSs 
ejected from the Galactic center. 

3.1. Background 

High Galactic latitude B-type stars were first reported 
by Humason & Zwicky (1947). Early spectroscopic stud- 
ies (Feige 1958; Berger 1963; Greenstein 1966) showed 
that high latitude B-type stars are a mix of Popula- 
tion II post-main sequence stars and Population I main 
sequence run-aways. Accurate space motions and dis- 
tances are now known for nearby d < 10 kpc B-type stars 
observed by Hipparcos. A detailed analysis by Martin 
(2004) shows that 2/3 of the high latitude B stars in the 
Hipparcos catalog are main sequence run-aways and 1 /3 
are evolved (mostly blue horizontal branch) halo stars. 
Martin (2004) finds that all of the main sequence run- 
away B stars in the Hipparcos catalog have orbits consis- 
tent with a disk origin. 

Our HVS survey probes much fainter and more dis- 
tant B-type stars than observed by Hipparcos. Thus 
we cannot rely on proper motions or parallaxes. Brown 
et al. (2006b) discuss the nature of the late B-type stars 
and conclude that the distribution of radial velocities 
and mctallicities suggests that they are likely a Galac- 
tic halo population of post-main sequence stars and/or 
blue stragglers. We now re-visit the possibility of main 
sequence run-aways in our sample. 

Theorists have shown that unbinding a stellar binary 
through supernova disruption (Blaauw 1961) or binary- 
binary interactions (Poveda et al. 1967) produces a maxi- 
mum ejection velocity of 200 - 300 km s _1 (Leonard 1991, 
1993; Portegies Zwart 2000; Davies et al. 2002; Gualan- 
dris et al. 2004; Dray et al. 2005). The maximum velocity 
is set by the escape velocity from the stellar surface. For 
example, a contact binary containing two 3 M Q stars has 
a Keplerian orbital velocity of 240 km s _1 . Extracting 
a large ejection velocity from disrupting such a system 
is difficult. Portegies Zwart (2000) considers supernovac 
unbinding binaries and finds that 90% of ejections have 
velocities between and 100 km s _1 ; only 1% of 3 M Q 
stars are ejected at 200 km s -1 . Davies et al. (2002) per- 
form similar calculations for binaries that undergo a com- 
mon envelope phase followed by a mass-transfer phase 




Fig. 3. — Galactic rest-frame velocity histogram of the clean 
sample of 560 B-type stars (upper panel). The best-fit Gaussian 
(thin line) has dispersion 105±5 km s . The lower panel plots the 
residuals of the observations from the best-fit Gaussian, normalized 
by the value of the Gaussian. In addition to the unbound HVSs, 
there is a significant asymmetry of fl positive velocity outliers 
+275 < v rf < 450 km s" 1 including SDSS J144955.58+310351.4. 



that produces a type la supernova. In such systems, 
Davies et al. find that 3 M Q secondaries have typical 
ejection velocities ^100 km s _1 and maximum ejection 
velocities of 250 km s _1 . By comparison, the escape ve- 
locity of the Galaxy near the Sun is at least 500 km s _1 
(Carney et al. 1988). 

Because main sequence run-aways travel on bound or- 
bits, they spend most of their time near the apex of their 
trajectories above the disk with small linc-of-sight veloci- 
ties. For example, a star ejected at 240 km s _1 vertically 
out of an infinite disk with surface mass density 95 M Q 
pc~ 2 reaches an apex of z = 11 kpc. This simple calcu- 
lation illustrates that main sequence run-aways cannot 
simultaneously have large velocities v > +275 km s _1 
and large distances z > 10 kpc. Davies et al. (2002) 
use a full Galaxy potential model and compute orbits 
for 1000 run-away stars randomly ejected from their pre- 
dicted velocity distribution. The vast majority of the 
Davies et al. (2002) run-aways are located at \z\ < 5 
kpc; at high Galactic latitudes \b\ > 30° the most dis- 
tant run-aways are found ^25 kpc from the disk. 

3.2. The Asymmetric Velocity Distribution 

Figure 3 plots the distribution of linc-of-sight veloci- 
ties, corrected to the Galactic rest-frame, for our clean 
sample of 560 late B-type stars plus HVS1 (Brown et al. 
2005a). It is apparent from Figure 3 that the velocities 
are well described by a Gaussian distribution. We itera- 
tively clip 3a outliers and calculate a 105 ±5 km s _1 dis- 
persion and a +14 ± 4 km s _1 mean for the distribution. 
We note that the full sample of 699 stars has a statis- 
tically identical distribution. The lower panel of Figure 
3 plots the residuals of the observations from this Gaus- 



4 



Brown et al. 



sian distribution, normalized by the value of the Gaus- 
sian. Stars with velocities \v r f\ < 275 km s _1 show low- 
significance deviations from a Gaussian distribution and 
probably constitute a halo population of post-main se- 
quence stars and/or blue stragglers (Brown et al. 2006b). 
The highest-significant outliers are the first HVS (HVS1), 
the four HVSs previously discovered in this survey (HVS4 
- HVS 7), and a new object: SDSS Jf 44955. 58+3f 0351.4. 

The possible HVS, SDSS Jf 44955. 58+3f 0351. 4, is a 
bright g' a = f 5.70 ± 0.03 star with an Af ±1 .4 spectral 
type and solar metallicity [Fe/H]^ = —0.3 ± 0.6. A 
main sequence star of this type has heliocentric distance 
d ~ 17 kpc (Schaller et al. 1992); an evolved blue hori- 
zontal branch (BHB) star has d ~ 7 kpc (Clewley et al. 
2005). Located at (l,b) = (48.7°, 63.9°), the star is high 
above the Galactic plane at z ~ 15 kpc and a Galacto- 
centric distance of R ~ 17 kpc (assuming d = 17 kpc). 
The star's +363 ± 10 km s" 1 heliocentric radial velocity 
corresponds to a minimum velocity in the Galactic rest 
frame of +447 km s _1 (sec Brown et al. 2006b). The star 
is listed with no proper motion in the USNOB1 (Monet 
et al. 2003) catalog, and thus is probably bound to the 
Milky Way. 

Curiously, we find a set of more distant stars with 
equally large, but bound, velocities. Figure 3 shows a 
significant asymmetry in the wings of the observed ve- 
locity distribution: If stars (excluding the 5 unbound 
HVSs) have large positive velocities +275 < v r f < +450 
km s _1 ; only 1 star has a comparable negative velocity 
at v r f = —286 ± 12 km s _1 . This asymmetry was first 
pointed out by Brown et al. (2006b) but is more signifi- 
cant in this larger sample. There is a 0.00022 probability 
of drawing 11 stars with +275 < v r f < 450 km s _1 and 
1 star with —450 < v r j < —275 km s _1 from a Gaussian 
distributions with the observed parameters. Integrating 
the wings of a Gaussian with the observed parameters, 
we would expect to find 3.6 stars with +275 < v r f < 450 
km s _1 and 1.6 stars -450 < v r j < -275 km s _1 . Thus 
there is an excess of ~7 stars with large positive veloc- 
ity and no apparent excess of stars with large negative 
velocity, significant at the ~3.5tr level. 

If the positive velocity outliers are main sequence stars, 
they are located at large distances 20 - 80 kpc (see Ta- 
ble 1). Yet the outliers cannot be run- away stars, be- 
cause the run-away mechanisms cannot produce simul- 
taneously large velocities and large distances. If, on the 
other hand, the positive velocity outliers are halo stars 
on radial orbits, we would expect to find equal numbers 
of stars moving toward and away from us, contrary to 
observation. Compact binary systems may also produce 
outliers in the velocity distribution. But velocity outliers 
resulting from compact binaries should be distributed 
symmetrically, again contrary to observation. We note 
that neutron stars are observed traveling at velocities of 
fOOO km s _1 , but such neutron stars are the remnants 
of asymmetric supernova explosions, and are not B-typc 
stars. 

3.3. Bound HVSs 

One explanation for the significant excess of B-type 
stars traveling +275 < v r f < 450 km s _1 is the HVS 
mechanism. HVSs are ejected by three-body interac- 
tions involving a MBH (Hills 1988). HVSs almost cer- 
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Fig. 4. — Predicted velocity distributions of 3-5 Mq HVSs 
with distance from the Galactic center (Bromley et al. 2006). We 
calculate distances for our B-type stars assuming they are main 
sequence stars: the four HVSs discovered by this survey (solid 
squares) and the 11 positive velocity outliers +275 < v r j < 450 
km s" 1 (open squares) including SDSS J144955. 58+310351.4 (open 
square with plus sign). 



tainly come from the Galactic Center. The presence of a 
3.6 x fO 6 M Q MBH (Ghez et al. 2005; Eisenhauer et al. 

2005) in the crowded Galactic center inevitably produces 
a "fountain" of HVSs ejections. The actual HVS ejection 
rate depends on the number of stars with radial orbits 
in the MBH's so-called "loss cone." If the MBH's loss 
cone is replenished with stars scattered in by massive 
star clusters, molecular clouds, or IMBHs (Perets et al. 

2006) , then the ejection rate can be orders of magnitude 
larger than the original Yu & Tremaine (2003) predic- 
tions. Regardless of the exact rate, this fountain picture 
results in a broad spectrum of HVS ejection velocities 
which include bound and unbound orbits (Hills 1991; Yu 
& Tremaine 2003; Gualandris et al. 2005; Levin 2006; 
Baumgardt et al. 2006; Perets et al. 2006; Sesana et al. 
2006; Bromley et al. 2006; O'Leary & Loeb 2006). Brom- 
ley et al. (2006) recently calculated ejection velocities for 
binaries containing 3-4 M Q primaries matched to our sur- 
vey of B-type stars. 

Over the volume sampled by our survey, Bromley et al. 
(2006) predict comparable numbers of HVSs ejected into 
bound and unbound orbits. Our survey has discovered 
4 HVSs on unbound orbits suggesting that 4 ± 2 of our 
excess ~7 positive velocity outliers are plausibly HVSs 
on bound orbits. Although the predicted and observed 
numbers of bound HVSs are statistically consistent, we 
note that additional HVS mechanisms may be at work. 
For example, O'Leary & Loeb (2006) show that single 
stars can be ejected by encounters with stellar mass black 
holes orbiting the central MBH. Such HVSs tend to have 
lower ejection velocities than HVSs ejected by the Hills 
(f988) mechanism, and thus may account for additional 
HVSs on bound orbits. 

Figure 4 plots the Bromley et al. (2006) predicted ve- 
locity distributions for main sequence HVSs ejected by 
the central MBH as a function of Galactocentric radius. 
HVS velocities increase with depth because only the most 
rapidly moving HVSs survive long enough to reach large 
Galactocentric distances. The large squares in Figure 
4 are our 4 HVSs (solid squares) and the f f stars with 
+275 < v r f < 450 km s -1 (open squares). We assume 
that the stars are main-sequence stars, and use the ob- 
served colors and spectral types of the stars to estimate 
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stellar luminosity, distance, and stellar mass (see Table 
1). The solid lines in Figure 4 indicate our survey limits 
for main sequence stars of a given stellar mass. 

We evaluate the probability of drawing our observed 
sample from the Bromley et al. (2006) predicted dis- 
tributions. We use the two-dimensional, two-sample 
Kolmogorov-Smirnov test (Press et al. 1992) to compare 
the observed and predicted distributions of HVS veloc- 
ity and distance. First, we generate 10 5 realizations of 
the observations - 4 HVSs plus 7 bound HVSs randomly 
drawn from the 11 high velocity outliers - and calculate 
the maximum difference in integrated probability against 
the Bromley et al. (2006) distributions. We use only the 
portion of the Bromley et al. (2006) distributions that 
fall within our observational survey limits (see Figure 
4). Second, we randomly sample 11 objects from the 
Bromley et al. (2006) distributions 10 5 times and calcu- 
late the maximum difference in integrated probabilities 
against the Bromley et al. (2006) distributions. Again, 
we use only the portion of the Bromley et al. (2006) dis- 
tributions that fall within our observational survey lim- 
its. Finally, the likelihood is the fraction of time that the 
probabilities of the synthetic data sets exceed the prob- 
abilities of the real data sets. The mean likelihood of 
drawing our observed sets of stars from the predicted dis- 
tributions is approximately 0.128, which neither strongly 
supports nor strongly rejects the Bromley et al. (2006) 
model. We conclude that the most plausible explanation 
for the observed excess of positive velocity outliers is the 
HVS mechanism ejecting stars from the Galactic center 
on bound orbits. 

We note that other stars in the literature may also be 
explained as bound HVSs. For example, the 5 M main 
sequence B star HIP 60350 has a heliocentric velocity of 
+230 km s _1 and a full space velocity of 417 km s _1 
based on Hipparcos measurements (Maitzen et al. 1998). 
The star's large radial motion U — —352 km means 
that it originated from well inside the solar circle, per- 
haps from the Galactic center. 

3.4. Mystery of the Early B-type Stars 

Our full sample of B-type stars contains many early- 
type stars at bright magnitudes. Four of the early B 
stars in the full sample are possibly bound HVSs with 
+275 < v r f < 400 km s _1 . If they are main sequence 
stars, the four early B stars have masses of 5 - 7 M Q 
and they are located at distances ranging 30 - 60 kpc 
(bounded by the FAST survey limiting magnitudes). We 
now ask whether the number of early and late B-type 
bound HVSs are consistent with a standard initial mass 
function. Assuming the stars are main sequence stars, 
there are nine 3-5 M stars in our sample located in 
the same volume as the four 5-7 M stars: a ratio of 
2.25. A Salpetcr initial mass function predicts a ratio of 
2.7, similar to what we observe. However, the lifetimes 
of 5 - 7 M Q main sequence stars are only 40 - 90 Myr 
(Schaller et al. 1992), a factor of 2 - 3 shorter than 3 - 
5 M stars. The travel times of the four early B stars, 
assuming v r f is their full space velocity, range from ^100 
to ~150 Myr from the Galactic center. Thus if the early 
B stars are main sequence stars ejected from the Galactic 
center, they do not have high enough velocities to survive 
to their inferred distances. 

How can we reconcile the lifetimes of the early B-type 
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a Galactocentric distance for a main sequence star of the observed 
spectral type. 

b Galactocentric distance for a BHB star of the observed color. 



stars with their travel times from the Galactic center? It 
is possible that the early B-type stars are blue stragglers, 
main sequence stars that have undergone mass transfer 
or mergers. Blue stragglers are usually associated with 
globular clusters, yet surveys of field BHB stars in the 
halo find that half of the stars are in fact high surface- 
gravity blue stragglers or A dwarfs (Norris & Hawkins 
1991; Preston et al. 1994; Wilhelm et al. 1999; Brown 
et al. 2003; Clewley et al. 2004; Brown et al. 2005b). 
Because HVSs likely originate in tight binary systems, 
it is possible that some HVSs experience mass-transfer 
from their former companion prior to being ejected from 
the Galactic center. A blue straggler origin is one expla- 
nation for HVS2, an 8 M HVS with a main sequence 
lifetime otherwise inconsistent with travel time from the 
Galactic center (Edelmann et al. 2005). 

It is also possible that the early B-type stars are main 
sequence stars ejected elsewhere in the Galaxy by dynam- 
ical interactions with compact objects like intermediate 
mass black holes (IMBHs). However, as the mass of a 
black hole decreases, the cross-section of HVS interaction 
decreases (e.g. Hills 1988). A stellar binary must come 
much closer to an IMBH before the gravitational tidal 
energy exceeds the binding energy of the binary. Thus 
an IMBH produces fewer HVSs than a MBH in the same 
environment. Unless IMBHs arc ubiquitous in dense stel- 
lar environments such as the centers of globular clusters, 
we expect that the overall HVS rate is dominated by the 
MBH in the Galactic center. 

A final possibility is that the early B-type stars are 
post main-sequence stars. Demarque & Virani (2007) 
argue that the Galactic center contains many more old 
dwarf stars than young B stars, and thus the MBH 
should eject many more BHB stars than main sequence 
B stars. Davies & King (2005) and Dray et al. (2006) 
argue that the Galactic center contains tidally-stripped 
evolved stars (in effect, BHB stars with small hydrogen 
envelopes) that mimic B stars. B-type BHB stars with 
hot effective temperatures and small hydrogen envelopes 
spend most of their helium-burning lifetimes at intrinsi- 
cally faint luminosities (Yi et al. 1997; Dray et al. 2006). 
Thus, if they are BHB stars, the early B-type stars are 
located 10 - 20 kpc from the Galactic center and have 
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Fig. 5. — Aitoff sky map, in Galactic coordinates, showing the clean sample of 560 B-type HVS candidates. Radial velocities, in the 
Galactic rest frame, are indicated by the color of the solid squares: purple is -300, green is 0, and red is +300 km s" 1 . The HVSs are 
completely off this color scale and are marked by: HVS1 (plus), the four HVSs discovered by this survey (solid squares), and the 11 positive 
velocity outliers +275 < v T j < 400 km s — 1 (open squares) including SDSS J144955. 58+310351. 4 (open square with plus sign). Dotted 
lines are lines of constant unbound HVS fraction calculated from the Bromley et al. (2006) models. 



travel times ranging from 25 - 50 Myr. The BHB pro- 
genitor lifetimes easily exceed the inferred travel times, 
reconciling any difficulties with the MBH ejection origin. 
To test the post-main sequence explanation requires high 
dispersion spectroscopy to determine the stellar nature 
of the bound HVSs. 

4. HVS SPATIAL ANISOTROPY 

The distribution of HVSs on the sky is interesting 
because the spatial distribution of HVSs is linked to 
their origin. Bromley et al. (2006) show that fraction 
of unbound HVSs increases with radial distance from the 
Galactic center (see Figure 4) . The fraction changes with 
depth because only the most rapidly moving HVSs sur- 
vive to reach large Galactocentric distances. Because the 
Sun is located 8 kpc from the Galactic center, our survey 
reaches approximately 16 kpc deeper towards the anti- 
center than towards the Galactic center. Thus, if the 
fountain model of HVSs ejected from the Galactic center 
is correct, we expect to find more unbound HVSs towards 
the anti-center. 

We use the Bromley et al. (2006) models (Figure 4) to 
calculate the expected fraction of unbound HVSs in our 
survey volume as a function of position on the sky. The 
predicted fraction of unbound HVSs ranges from 53% 
towards the Galactic center to 63% towards the anti- 
center. The fraction does not vary much because it is 
dominated by the survey volume at large distances. The 
HVSs found in this survey, however, are not located at 
distances beyond ~80 kpc. If we limit the survey vol- 
ume to the observed range of HVS distances, then the 
predicted fraction of unbound HVSs ranges from 31% to 
46%. We now consider the observations. 

Figure 5 plots the distribution of our clean sample 
of 560 B-type stars across the sky in Galactic coordi- 
nates. Color indicates the radial velocity, in the Galac- 



tic rest frame, of the stars. The four unbound HVSs 
(solid black squares) discovered by this survey are com- 
pletely off the color scale; HVS1 is marked by a plus 
sign. The 11 possible bound HVSs are indicated by open 
squares; SDSS J144955. 58+310351.4 is plotted with an 
open square with a plus sign. The dotted lines in Fig- 
ure 5 are lines of constant HVS fraction calculated from 
the Bromley et al. (2006) models. We note that our sur- 
vey is currently incomplete in the SDSS equatorial slices, 
located at I > 210° and b < 0°. 

It is striking that all of the unbound HVSs in our sur- 
vey are located towards the anti-center. We randomize 
the observed velocities among the positions of all 560 
stars, and find a 0.064 probability of randomly draw- 
ing the 4 unbound HVSs from our survey in the hemi- 
sphere 90 < I < 270°. Thus the unbound HVSs have a 
marginally anisotropic spatial distribution, being prefer- 
entially located in the anti-center hemisphere with ^2a 
confidence. We note that the 3 unbound HVSs not from 
this survey, HVS1 (Brown et al. 2005a), HVS 2 (Edel- 
mann et al. 2005), and HVS3 (Hirsch et al. 2005), are 
also located in the anti-center hemisphere 90 < I < 270°. 

The 11 possible bound HVSs, on the other hand, ap- 
pear more evenly distributed across the sky: 6 bound 
HVSs are located 90 < I < 270° and 5 bound HVSs 
are located \l\ < 90°. To calculate the variation in un- 
bound HVS fraction, we consider only 7 bound HVSs 
(the excess number of high velocity outliers in Figure 3) 
and split them 4/3 between the anti-center and Galactic- 
center hemispheres, respectively. The fraction of un- 
bound HVSs is then 50±30% in the anti-center hemi- 
sphere and 0±40% in the Galactic-center hemisphere. 
Although small number statistics overwhelm our mea- 
surement, the unbound HVS fraction appears consistent 
with the predictions of Bromley et al. (2006) and the 
Galactic center ejection picture. 
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There is another possible explanation for spatial 
anisotropy: HVSs ejected via single star encounters 
with a MBH binary are preferentially ejected in the 
orbital plane of the MBH binary (Gualandris et al. 
2005; Holley-Bockelmann et al. 2005; Levin 2006; Sesana 
et al. 2006; Merritt 2006; Baumgardt et al. 2006). Re- 
cent N-body simulations for the case of an IMBH inspi- 
ralling into a MBH, however, show that HVSs are ejected 
nearly isotropically because the angular momentum vec- 
tor of the IMBH rapidly changes ona~l Myr timescale 
(Baumgardt et al. 2006). Moreover, our observed set of 
HVSs do not share common travel times from the Galac- 
tic center. Different IMBH inspiral events would likely 
have different orbital histories, thus it seems unlikely that 
the observed HVS spatial distribution can point back to a 
binary MBH origin. A set of HVSs with common travel 
times must be found in order to test the binary MBH 
origin. 

5. CONCLUSIONS 

We discuss our targeted survey for HVSs, a spectro- 
scopic survey of stars with B-type colors. We introduce 
a new component of the survey that samples brighter and 
bluer stars using the FAST spectrograph. Our combined 
set of observations now contains 699 stars of B spectral 
type, 188 white dwarfs, and 1 previously unidentified 
z = 2.43 quasar located at 15 h 58 m 5l!85, +22 d 21 m 59!9 
(J2000). 

We find no new unbound HVS, but we find evidence 
for a new class of bound HVSs. The velocity distribution 
of the B-type stars shows a significant excess of ~7 stars 
with large positive velocities +275 < v r f < 450 km s _1 
and no apparent excess of stars with large negative veloc- 
ity, significant at the ~3.5ct level. Neither halo stars on 
radial orbits, compact binary stars, nor main sequence 
run-away stars ejected from the disk can explain simul- 
taneously the large velocities and distances of these stars. 
The most plausible explanation is that the positive ve- 
locity outliers are HVSs ejected from the Galactic center 
on bound orbits. 

Our bright FAST survey contains a large number of 
early B-type stars, four of which have velocities +275 < 
v r f < 450 km s . These four possible bound HVSs have 
main sequence lifetimes too short to survive to their in- 
ferred distances, however. One explanation is that the 
stars are post-main sequence stars. Establishing HVSs 
as main sequence or post-main sequence stars is impor- 
tant for measuring the stellar mass function of HVSs 
and probing the types of stars that orbit near the MBH. 
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Stellar rotation is a useful discriminant between rapidly 
rotating main sequence B stars (Abt et al. 2002; Mar- 
tin 2004) and slowly rotating BHB stars (Peterson et al. 
1995; Behr 2003). Metallicity is also a good discriminant 
between young main sequence B stars and post-main se- 
quence halo stars. Echelle observations of the brightest 
HVSs, underway now, will reveal the stars' true nature. 

A Galactic center origin predicts that unbound HVSs 
should be found preferentially towards the anti-center as 
a result of survey selection. We observe that the unbound 
HVSs have a marginally anisotropic distribution at 2- 
o confidence. Our estimates of unbound HVS fraction 
are consistent with the Galactic center ejection picture. 
Measuring the spatial distribution of a larger sample of 
HVSs may provide a strong test of the stars' origin. 
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cies by registering background galaxies across Advanced 
Camera for Surveys images. The HVSs have predicted 
proper motions of 0.5 - 1 mas yr _1 , thus 3cr proper mo- 
tion measurements are possible over a three year baseline 
(Gnedin et al. 2005). 
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interactions with the central MBH. 
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APPENDIX 
DATA TABLE 

Table 2 [on-line data table, attached as tab2.dat in the source] lists the 643 new observations from our HVS survey, 
excluding the extra-galactic objects. 134 HVS candidates were observed with the MMT and 509 were observed with 
FAST (see §2). Table 2 includes columns of RA and Dec coordinates (J2000), g' apparent magnitude, (u' — g') and 
(g 1 — r')o color, and our heliocentric velocity Vheiio- 

The column WD indicates whether the object is a B-type star 

(WD=0) or a white dwarf (WD=1). 

We note that (u' — g')o color correlates strongly with the observed spectral type of our B-type stars. Our spectral 
types are completely independent of the SDSS photometry. We classify spectral types based on O'Connell (1973) and 
Worthey et al. (1994) line indices as described in Brown ct al. (2003). Typical uncertainty is ±1.6 spectral sub-types. 
Figure A6 plots our spectral types and the (u' — g')o color for our full sample of 699 B-type stars. The solid line is the 
best-fit relation: spectral type = 9.5(u' — g') + 11.5, where spectral type 10=B0, 15=B5, 20=A0, etc. The spectral 
types of the stars in our HVS survey range B3 - A2, with a dispersion around the best-fit relation of ±1 spectral 
sub-types or ±0.1 mag in (u' — g') . 
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Fig. A6. — Spectral types determined from our spectroscopy 
and SDSS («' — g')o color for the full sample of 699 B-type stars. 



